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Abstract 

Molecular motors are single macromolecules that generate forces at the pi- 
conewton range and nanometer scale. They convert chemical energy into me- 
chanical work by moving along filamentous structures. In this paper, we study 
the velocity of two-head molecular motors in the framework of a mechanochem- 
ical network theory. The network model, a generalization of the recently work 
of Liepelt and Lipowsky (PRL 98, 258102 (2007)), is based on the discrete 
mechanochemical states of a molecular motor with multiple cycles. By general- 
izing the mathematical method developed by Fisher and Kolomeisky for single 
cycle motor (PNAS(2001) 98(14) P7748-7753), we are able to obtain an exphcit 
formula for the velocity of a molecular motor. 
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1 Introduction 



In biological cells, molecular motors are individual protein molecules that are 
responsible for many of the biophysical functions of the cellular movement and me- 
chanics. Important examples of motor proteins are kinesin [H [21 E], dynein [U |3], 
mysion O [71 E] and Fq-Fi-ATP synthase [9j. Molecular motors are mechanochemical 
force generators which convert biochemical energy (stored as ATP, adenosine triphos- 
phate) into mechanical work in a thermal environment [IDl [H]. Many molecular 
motors, due to their two-head nature and hand-over-hand mechanism, can move pro- 
cessively along their tracks for a long time before its dissociation from the track. For 
example, myosin slides along an actin filament, kinesin and dynein along microtubule 
(MT). The velocity of molecular motors is quite fast, with mean velocity at about sev- 
eral hundreds nanometers per second [12] . Understanding how the various molecular 
motors operate is a significant scientific challenge with important nano-engineering 
implications. 

To understand the principle of molecular motors, a good mathematical model is 
essential. Much progress has been made in recent years in theoretical analysis of 
molecular motors. Mainly two different approaches have been taken: The ratchet 
models that consider motor chemical transitions occur without explicit coupling to 
motor steppings [l3l[Tlj, and the discrete chemical models that contain only a single 
chemomechanical cycle [151 HB]- Recently, however, Liepelt and Lipowsky [T71 [T8] 
introduced a six-state network to model the chemomechanical motor cycles, in which 
the dynamics of two-head motor molecule is described by a Markovian jump process. 
In [19], Schmiedl and Seifert used a two states network to discuss the efficiency of 
the molecular motors. The importance of the latter development is in introducing 
futile cycles into the discrete chemical model, thus making the discrete chemical 
approach and continuous Brownian approach more connected. Their results indicate 
that the network modeling approach is a good choice for the theoretical analysis of the 
molecular motors. In the past, a great deal of mathematical analysis is based on the 
Brownian ratchet formalism. Similar network models has also be used successfully in 
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the theoretical analysis of other biochemical processes [20l [2T] . 

In this paper, we shall generalize the network model to include arbitrary 2N num- 
ber of states. In particular, we shall use the network model to analyze the movement 
of molecular motors. Mathematically, therefore, the models developed in [iHl [19] and 
even those in [201 El] can be regarded as special cases of our network model. In the 
framework of this network model, we further develop a method pioneered by Der- 
rida. Fisher and Kolomeisky [22l [23l [23] to compute the mean velocity of a molecular 
motor. 

In our model, a two-head molecular motor with hand-over-hand mechanism is as- 
sumed to have 2N mechanochemical states in their movement, denoted by 0, 1, 2, ■ ■ ■ , 2N— 
1 (see Figure [1]). The two heads moves exactly with half of cycle out of phase. If 
there are 2N states in the hydrolysis kinetic cycle of a single head; we have states 
(0, A^), (1, + 1), (2, + 2), ■ ■ ■ , (2A^ - 2, - 2), and (2A^ - 1, A^ - 1) for the motor 
with two heads. The hand-over-hand mechanism means the motor "walks" a step 
with the transition (A^, 0) (0,A^), switching the leading and the trailing head. 
However, it is possible that the translocation does not occur, and the kinetic cycle is 
completed as a futile cycle, with two ATP hydrolyzed, one by each head. 

From now one, we shall use the state of the leading head to denote the state of the 
motor; and denote the forward and backward rate parameters at state i as Ui (i.e., 
i ^ i + 1) and Wi (i — > i — 1) respectively, which satisfy U2N+i = Ui and W2N+i = Wi 
(since molecular motors move forward periodically). Generally, the transition rates 
Ui and Wi depend on the external force F^xt and the free energy AG released by the 
fuel molecular. The transition rates between state A^ and 0, the hand-over-hand, 
are u and w. In the following, we suppose that all these transition rates are known 
explicitly. 

The transition from A^ — represents the switching between the leading and 
trailing heads, thus moves one motor step. If a mechanochemical process takes 
1^---^A^— s>A^+l— s> 2A^ — 1^0, the molecular motors make no 
mechanical step while hydrolyzing two ATP. However, if the process takes 0^1^ 
■ >^A^— >0— >---A^— >0, then the motor hydrolyzed two ATP and moved 

3 




Figure 1: A schematic depiction of the 2N states network model for molecular motors. 
One forward step of molecular motors is completed only in the biochemical process 

— i> 0. In mechanochemical process O^l^---— s>A^— s>A^ + l^---— > 2N — 1 
0, the molecular motors make no mechanical step while hydrolyzing two ATP. 

two steps. It can be readily found that, for A^ = 3, this model reduces to the 6 states 
network model in [T7|, for A^ = 1, this model reduces to the 2 states model in P^ . 

In the next section, we shall give the formulation of the velocity of molecular 
motors using the network model. We will discuss some special cases in section 3. The 
force dependence of the transition rates Ui,Wi and u^w is discussed in section 4. In 
section 5, we will discuss the continuous mechanochemical sate case of our multi-cycle 
model, and section 6 contains concluding remarks. 



2 The velocity of molecular motors 

In this section, we will calculate the velocity of the molecular motors in the frame- 
work of our network model. The method used in the following is similar to the one 
used in [221 [23] and [21] . 

Let pi{t) be the probability density for finding molecular motors in state i at time 
t. The evolution of the probability density pi{t) is governed by the following master 
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equations 



dpi 
dt 



--{pi^lUi^l + Pi+iWi+i) - Pi{Ui + Wi) 
--{pi-lUi_i - PiWi) - {piUi - Pi+iWi+i) 



and 



= Ji - Ji+i for 1 < z < AT - 1 or iV + 1 < z < 2Ar - 1 



^ =iP2N-lU2N-l - PqWo) - (PqUq - PlWi) + {pnU - Pqw) 



dp 



•N 



-J2N — Jl + J 



--{Pn-iUn-i - PnWn) - {pnUn - Pn+iWn+i) - {pnu - Pow) 



dt 

—Jn — Jn+1 — J 

where 

Ji = pi-iUi-i - piWi J = PnU - pow 

Ji is the probabihty flux from mechanochemical state i — 1 to state i, and J 
probabihty flux from mechanochemical state N to state 0. At steady state, 



Ji = J2 



Jn Jn+1 — J. 



N+2 



J2N Jl — J2N + J 



By Eqs. dlH), one can know that 



i=l 

PN+k = PN 



Uj-l 



fe-1 k-l 



i=l j=i 



n 



i=l 



UN+i-1 
WN+i 



where 1 < A; < — 1 and 



Pn = Po 



N 

n 

i=l 



Uj-l 
Wi 



1=1 j=i 



N-l N~l 

1=1 j=t 



Wk 

UN+j 
WN+j 



Wi 



Jl-J 

WN+k 



A 

Wn 



J2N = Jl — J means P2N-1U2N-1 — PqWq = Ji — J, so 
P2N-1U2N-1 Jl — J 



Po 



--Pn 



Wo 

N 

n 

i=l 



Wq 



Un+i-I 
WN+i 



N~1N~1 
1 = 1 j=x 



UN+j 
WN+j 



Jl - (pnu - pow) 



Wo 



Substituting into ([7j), one obtains 



where 



A 



and 



2N-1 2N-1 



i=l j=i ^ ■' 



1 

— + 



Po = ^Ji 



Af-l N-1 



B 



Wn 



1=0 



U T-i- [ Uj \ w 
m^T A J. \ 7/1. / 



i=l j=i 



2N-1 2N-1 



2N-1 2N-1 



i=N+l j=i ^ -' 



U 



^-5: n ^7 

i=Ar+l j,-=i 



n 

i=0 



— 1-1 



WqWn 

(9) 



(10) 



So 



J 



N 



U 



n 



Wi 



w 



A 
B 



N-l N-l 



U 



— )Ji= : CJi 



(11) 



By dS]) (jS]) dHD (fTTl) . we get the expressions of probabities pk and Pat+j; as functions of 



Pk 



PN+k 



Awt 
B 

Awq 
B 



k-l 



i=Q 
N+k-1 

n 



k-l k-l 



i=l j=i 



Jl 



i=0 

N+k~l N+k-1 



N+k-1 N+k-1 



s n (2 

1=1 j=i 



n 

i=N+l j=i 

Conservation of probabihty requires 

27V-1 

5^ Pfc = 1 



(1-C) 



Jl 



WN+k 



k=0 



So, from (fll) (til ((11), one knows 



(12) 



(13) 



(14) 



Jl 



D 



where 



A f'i^' 



k=l 
2N-1 



k-l 

-n 

?.=0 



k=l 



Wi 
2N-1 



2N-1 



1 -E 



k-l k-l 



J-1 \Wk) ."^.Wk 



k=N+l 



k=l 
k-l k-l 



Wk ^ \Wj 

1=1 j=i ^ ■' 



i=N+l j=i ^ ■' 



(15) 



(16) 



In view of ( ITTj) and (fT5|) . one obtains 



J = CJ, = ^ (17) 



So the mean velocity of the molecular motors is 

CL 



where L is the stepsize of the molecular motors (8.2nm for motor protein kinesin). 
Certainly, the explicit expresions of probabilities Pk{0 k < 2N — 1) also can be 
obtained bv (|6|) (|8l) (fT2l) (fT3l). 



3 The special cases of the network model 



In this section, we consider some special cases of the network model. 
(1) Wq = W2N-1 = W2N-2 = ■ ■ ■ = wn+1 = 0, and un = mat+i = ■ ■ ■ = U2N-1 = (see 
Figured (Up)): 

In this case, our network model reduces to the usual one dimensional hopping 
model [211 [23] • It can be easily found that Jj = pi^iUi^i— piWi = for < i < 2N, 



and J = pnu — PqW = PqUq — piWi = Ji- By Ji = J2 



pk = po 



n 

i=l 



Uj-l 
Wi 



k-l fc-1 



■ 1 ■ ■ x^i 



Jni one obtains 
J_ 



(19) 



and pn = Pollti (^) - [1 + Eti'nf=/ (^)] ^- At the same time, p^u 
PqW = J implies 

N / ^ r N-lN-l 

PN = PoYl 



J + PqW TT / Ui-1 



U 



i=l 



Wi 



J_ 

Wn 



— (20) 



which gives 



Po 



u 



nN ( Ui-i 
1 = 1 \ Wi 



J 



(21) 



w 



Combing (HHD ([HD, we get 

T-riV-l / u. 



1 



Pk 



1 + e:li^ 



J=t \w 



nN ( Ui^i 
t = l \ Wi 



W 



n 

i=l 



Wi 



fe-1 k-1 

1=1 j=i 



Wk 



J 



(22) 




Figure 2: Special cases of the network model: (Up) in which Wq = W2N-1 = W2N-2 = 
■ ■ ■ = Wn+1 = 0, and Un = u^+i = ■ ■ ■ = U2N-1 = 0. (Down) in which wq = Wi = 
■■■ = W2N-1 = 0. 
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Finally, Ylk=o Pk = ^ gives 



where 



A 



1 + 






u 




"n.1:( 


^-^ — W 
m J 





J 



N 

E 

k=0 



(23) 



n 

i=l 



Uj-l 
Wi 



N 

E 

k=0 



fc-1 A;-l 



1=1 j=« ^ - 



1 



(24) 



So in this case, the mean velocity of molecular motors is = JL = L/A, and the 

probabilities pk are given by Eqs. fl2T]) fl22|) . 

(2) = Wi = ■ ■ • = W2Af~i = (see Figure [2] (Down)): 

In this case, Jj = pi^iUi-i for 1 < i < 2A^, and J = pnu — pow. At the steady 
state 

Pk = —Po PN+k = -^^Pn for < A; < - 1 (25) 

Uk UN+k 

Due to Jn+i + J= J + N — \, one knows 

Pn{u + un) = Po{uq + w) (26) 



I.e. 



Uq + W 

PN = ■ Po 

U + Un 



in view of fl25|) fl27|) and J2k=o Pk = 1; obtains 

1 



Po 



hence 



'^(Mo + t^) 
J =pArM - poW = ■ W ] Po 

U + Un 

UUq — U nW 



Mo(m + un) X;f=o^ i + un{uo + w)Y1 



2N~1 J_ 
^k=N Uk 



and the probabilities pk can be obtained by Eqs. ( |25l) (l27j) ( l28l) . 
(3) = 1: 
In this case, 



(27) 



(28) 



(29) 



A = r> = 



(30) 
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f UUo \ f Ui _i U \ f UUq W UqUi \ u 

C = w ] { hwo H \ 1 

\WiWo WiWoJ \Wi ) \WiWo Wo WiWo ) 

and 

UUQ JiL+ "0"! -1 WN 

w-\ wn wi) w-\ wr 



The probability fliix of the special case (1) is 



J = ^f!fo^^f!!f!^ (32) 

Uo ^ W\ ^ u ^ w 



The probability flux of the special case (2) is 



J = (33) 

U + Ui + Mo + W 



(4) AT = 2: 

In this case, the probability flux ([23D is 



^ WiW2W'iWo 



+ U2U?^W\ + W1W2W3 + W1W2US + UW1W3 + UW1U3 + UU1W3 + UU1U3 

(34) 



the probability flux (l29l) is 

J _ juUo - U2W) U1U3 ^g^^ 

UU1U3 + U3UU0 + U1U2U3 + U3U0U2 + U1U0U3 + U1U0U2 + U1WU3 + U1U2W 



4 The force dependence of the transition rates 

As pointed out in the introduction, the transition rates Wj, Wj, m, w depend on the 
external force F . For nonzero external force F, the force dependence of the transition 
rates Ui,Wi,u,w can be modeled as the following 

(36) 

Ui = fc+e"^^"^^' Wi+i = kr^y^^'^^-^^^^ < z < 2iV - 1 
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chemical state 



N H+S 



Clteiuical stale 



Figure 3: Energy profile of a molecular motor in the neighborhood of local equilib- 
rium mechanochemical state: (Left) Molecular motor undergoes thermal fluctuations 
around the i— th local equilibrium position with potential Ei, which corresponds to 
mechanochemical state i. It moves forward (to the right) or backward (to the left) 
when it acquires enough energy to across the energy barriers Ej^ or Ef_^. The local 
equilibrium position i and i + 1 are separated by characteristic distance Lj, the local 
equilibrium state i and the transition state i + 5i with energy Ef are separated by 
characteristic distance SiLi, and the local equilibrium state i + 1 and the transition 
state i + 5i are separated by characteristic distance {1 — Si) Li. The energy difference 
between state i and i + 1 is AEi = Ei — Ei+i. (Right) Molecular motor undergoes 
thermal fluctuations around the A'"— th local equilibrium position with potential E^, 
which corresponds to mechanochemical state N. It moves forward (to the right) when 
it acquires enough energy to across the energy barriers Ej. The local equilibrium N 
and are separated by characteristic distance Ls, the local equilibrium N and the 
transition state N + S with energy Ej are separated by characteristic distance 5Ls, 
the local equilibrium and the transition state N + 6 are separated by characteristic 
distance {1 — 5)Ls. The energy difference between state and N is A// — AE. 
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where j3 = l/ksT, < S,Si < 1 are load distribution factors that reflect how the 
external force affects the individual rates (THl [221 121] (see Figure [3]), Lq + Li + ■ ■ ■ + 
Ln-1 + Ls = L, Lq + Li + ■ ■ ■ + Ln^i = Ln + Ln+i + ■ ■ ■ + L2n-i- 

In (!36l) . the load distribution factors 5 and 5i can be determined by experimental 
data as in [251 EEl [271 [2H] • Thermodynamic consistency requires kf / k~^^ = e^^^* and 
k^ /k^ = e^^^^^^^\ where AEi = Ei—Ei^i is the potential energy difference between 
mechanochemical states i and i + 1 (see Figure [3D, AE = J2il=o^ = Eq — E^ 
is the potential energy difference between mechanochemical states and state A^, in 
the no external force case, which is the energy barrier of the movement of molecular 
motors. A/i is the chemical energy transferred to the motors in one mechachemical 
step, which comes from the hydrolysis of the fuel molecule ATP (see Figure [1]). 

5 Continuous mechanochemical state multi-cycle 
network model 

As the number of mechanochemical states 2A^ tends to infinite, our multi-cycle 
network model (see Figure [ID can be approximated by the continuous mechanochemi- 
cal state model (see Figure[5l). In this model, there're two chemical passway from state 
Ofe to state Ik with different potentials Vi(x) and V2{x) {kL < x < kL+a) respectively. 
From state 1^ to state 0^+1, the potential function is V^^x) {kL + a < x < {k + l)L). 
Biophysically, the potentials Vi{x) are periodical, i.e Vi{x + L) = V{x), and satisfy 
^i(Ofc) = V2{0k),V^{U) = V^iU) = V,{h). 

In the i—th chemical passway, the motion of molecular motors can be described 

by the following Fokker-Planck equation: 

Pi{x,t) _ Piix,t) dVi{x) ^ j^ dpi{x,t) \ _ dJj 
dt (9x \ 7 dx dx J dx 

where kL < x < kL + a for i = 1,2 

(37) 

kL + a < X < {k + 1)L for i = 3 
in which 7 is viscous friction coefficient, D is free diffusion coefficient which satisfies 
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Chemical state 



0(2N3 N 0(2N) 



Figure 4: The energy profile in mechanochemical cycles: (Left) No external force 
case: the energy barrier between mechanochemical states and N is |A£'|. Af- 
ter mechanochemical state A^, the molecular motor might back to state through 
mechanochemical passway A^^A^ — l^---^OorA^^A^-|-l^---^ 2A^(0). In 
this case, the molecular motor makes no any mechanical steps. Also, the molecular 
motor might back to state immediately through mechanochemical passway A^ — 0. 
In such case, molecular motor completes one forward mechanical step, with one fuel 
molecule ATP is hydrolyzed. The free energy released by one ATP molecule is A//. 
(Right) Nonzero external force F case: in this case, the energy barrier between 
mechanochemical states and A" is | A£J| + F{L — L^), which is bigger than the no 
external force case. So it will be more difficult for molecular motors to make a forward 
step. During one forward step, the energy dissipation is A/x — FL, which is small 
than the no external force case, since part of the energy A// released by the ATP is 
used to do useful mechanical work. 
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Figure 5: Depiction of continuous mechanochemical state multi-cycle network model: 
There're two chemical passway between mechanochemical state Ofc and Ifc, in which 
the potentials are Vi{x) and V2{x) {kL < x < kL + a) respectively. The potential 
between mechanochemical state and Ofe+i is Vslx) {kL + a < x < {k + 1)L). 

Einstein relation D — ksT/^, Pi{x,t) is probability density for finding molecular 
motors at mechanochemical state x in i—th passway at time t and Ji{x,t) is the 
probability flux. Deflne 

oo oo 

p^{x,t) = p^{x + kL,t) Ji{x,t) = Ji{x + kL,t) (38) 

fe=— oo fe=— oo 

it can be readily verified that 

pi{x,t) ^ Pt{x,t) dVi{x) _^ ^ dpi{x,t) \ ^ dJj 
dt dx \ 'y dx dx J dx 

where < x < a for i — 1,2 

(39) 

a < X < L for i = 3 

At steady state, the probability fiux Jj is constant and the probability Pi{x) satisfies 

dpi{x,t) _^ dVi{x) Pi{x,t) _ _Ji, 
dx dx ksT D 

where 0<x <a for i = 1, 2 (40) 
a < X < L for i — 3 
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Under the following constraints 



i pi{x)doi 
Jo 



x+ p2{x)dx + / p3{x)dx = 1 



'0 (41) 
pi{0) = p2{0) = psiL) pi(a) = p2(a) = P3(a) Ji + J2 = J3 



we can get the solutions of fj40|) 



D 

J. 



c^bt dy \ e ''bt for i = 1,2 



P3ix)=[Cs-^ I e^B 



(42) 



''"^ dy j e ^bt 
where the constants Ci, C2, C3 are the following 



Ci — C2 



/o f e '^s^ + e ^bt \ dy [ e^BTdy\ + ( e "b^ dy j ij^e ^bt dy j 



V3(L)-Vl(0) 

e ''bt a 



„a / '^it^) y2(i/) \ 1 / ^3(v) \ Vi(0)-V3{L) / V"i(;/) \ / pa ^2(1/) 

Jq I e ''bT + e ''sT' I (iy I e '=sr j g ^sT _|_ I g fe^T j I e ''sr dy 



A 



and the probability fluxes Ji, J2, J3 are 

/ Vi(0)-V3(£) 

Ji = I e - 1 

V'2(0)-'^3(^) 

J2 ■ 



' V2(l/) 

e^sT dy] D / A 
e '=B'r c/y D / A 



(43) 



■h = Ji + J2 
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The expression A is 
A = 



Vi(y) 

e dy 



°- V2(y) 



+ 

X 

+ 
+ 
+ 



'dy 



e "B 

L 



dy 



e ''bt dy 



e fe_gT 



dy 



e 



dy 



Viiy) 



V2iy) 



e ''bt dy + e '^s^ dy + e '^s^ dy 



V3(y) 



vi{y) 
e ^bt dy 



■ V2(y) 

e ''bt dy 



' V3(y) 
e ^bt dy 




3 fcpjT 



L ex 



'0 ^0 
fa i-x 




dydx 



Vi(y)-vi(.^) 

e ^bt dydx + 




dydx 



I a J a 



JO 




V-l(0)-V3(L) \ 

X ( 1 - e "bt 1 



Vliv) 



V2(v) 



+ 
+ 
+ 
+ 
+ 

X 



e ''bt dy ] I e dy 



• V2(v)~V2{x) 

e ^bt dydx + 

Viiy)-Vi(^) 

e '^bt dydx + 



Vs(y)-v^i(^) 

e '^bt dydx 



V3(y) 



e ''bt dy 



■ V2(y) 
e ''bt dy 




e B 

L rx 



'0 JQ 

fa nx 




''bT dy 

V3(i;)-V3(^) 

fcRT 



dydx 




'a J a 
-L rx 



JO 




Vi(y)-v-i(=^) 

e '^bt dydx 




V2{y)-V2{^) 

e ^bt dydx + 



L rL 



Jx 

■a ro, 




V3[y)-V3{^) 



dydx 




Vi{y)-Vi{^) 

e ^bt dydx + 



'a J x 



J X 
a 




V3(y)-V3(=') 



'dydx 



V2(y) 
e ^bt dy 

Vi(y) 



V3(y) 



e'^BT dy ] + [ e ''bt dy 



e ''bt dy 



• V2iy) 

e ''bt dy 



vi{y) 

e ^bt dy 



V2(y) 

e ^bt dy 



e ^bt dy 



e ''bt dy -\- e ''bt dy 



) kgT 



>0 



Therefore, in the framework of this continuous mechanochemical state multi-cycle 
network model, the expression of the mean velocity of molecular motors is 



V = {Ji + J2)a + J^{L - a) 



V-l(0)-V^;(£) 



/o (e'^BT + e '^bt J dy 



DL 



A 



(44) 



Obviously, V > if l^i(O) > V^^L) and < if ^1(0) < V^{L). It can be readily 
verified that the equations ^ ^ can be obtained by applying spatial discretization 
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to Fokker-Planck equation (1371) . with some detailed expression of the transition rate 
Ui, Wi (see 

6 Concluding remarks 

In this paper, a general multi-cycle network model of molecular motors is theo- 
retically discussed. The explicit formulation of the velocity has been obtained. This 
model can be regarded as a generalization of the one designed by Liepelt and Lipowsky 
in [TS] and the one used by T. Schmiedl and U. Seifert in [IS]. The method used in this 
paper is similar as the methods used by Derrida, Fisher and Kolomeisky [22 | [23l [2^ . 
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